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Abstract. We consider a two-dimensional geometrically frustrated integer-spin Heisenberg system that 
admits an exact ground state. The system corresponds to a decorated square lattice with two coupling 
constants Ji and J2, and it can be understood as a generalized Shastry-Sutherland model. Main elements 

O 

of the spin model are suitably coupled antiferromagnetic spin trimers with integer spin quantum numbers 
s and their ground state # will be the product state of the local singlet ground states of the trimers. 
We provide exact numerical data for finite lattices as well as analytical considerations to estimate the 

m 

range of the existence in dependence on the ratio of the two couplings constants J2 and Ji and on the 

■ spin quantum number s. Moreover, we find that the magnetization curves as a function of the applied 

o: 

£N) . magnetic field shows plateaus and jumps. In the classical limit s — > 00 the model exhibits phases of three- 

and two-dimensional ground states separated by a one-dimensional (collinear) plateau state at 1/3 of the 

■ saturation magnetization. 

a; 

1 Introduction spin systems is a challenging task. Exact statements about 

the properties of quantum spin system are known only in 



The concept of frustration plays an important role in the 
search for novel quantum states of condensed matter, see, 
e.g., [1I2I3I4I5] . The investigation of frustrating quantum 



exceptional cases. The simplest known exact eigenstate is 
the fully polarized ferromagnetic state. Furthermore the 

~ ; \ one- and two-magnon excitations above the fully polar- 
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see, e.g., [6171819] . An example for non-trivial eigenstates fragments of finite length can lead to an infinite series of 

is Bcthe's famous solution for the one-dimensional (ID) magnetization plateaus [22] . 

Heiscnbcrg antiferromagnet (HAFM) [TO]. The investiga- 0thcr examples of product ground states are singlc- 

tion of strongly frustrated magnetic systems surprisingly gpin product statcs f 1D X YZ model [M] and the the 

led to the discovery of several new exact eigenstates. Some highly deg enerate ground-state manifold of localized-magnon 

of the eigenstates found for frustrated quantum magnets states found for ant if e rromagnetic quantum spin systems 

arc of quite simple nature and for several physical quanti- on various frustrated lattices [25]. Finally, we mention the 

ties, e.g., the spin correlation functions, analytical expres- so _ ca n ed central-spin model or Heisenberg star where also 

sions can be found. Hence such exact eigenstates may play cxact sta t cmC nts on the groundstate are known gg]. 
an important role either as groundstates of real quantum 

Although, at first glance such singlet-product states 

magnets or at least as groundstates of idealized models . . . 

seem to exist only for exotic lattice models, it turned 

which can be used as reference states for more complex 

out that such models are not only a playground of thc- 

quantum spin systems. A well-known class of exact eigen- 

oreticians but may become relevant for experimental re- 
states are dimerized singlet states, where a direct product . . 

search. The most prominent example is the above men- 

of pair singlet states is an cigenstate of the quantum spin __ 

tioncd Shastry-Sutherland model introduced m 1981 [12] 

system. Such states become groundstates for certain val- 

for which only in 1999 the corresponding quasi-two-dimen- 

ucs/rcgions of frustration. The most prominent examples 

sional compound SrCu2(B03)2 was found [27128] . Other 

are the Majumdar-Gosh state of the ID Ji — J2 spin- , . , 

examples are the quasi-lD spm-Peicrls compound CuGeOs, 

half HAFM [TT] and the orthogonal dimer state of the . . 

see, e.g., [29], or the star-lattice compound [Fc3(fi3-0)[fiT- 

Shastry-Sutherland model, see, e.g., |12|13|14|15|16|17|. _____ 

OAc) 6 -(H20)3][Fe3(^3-0)(Ai-OAc) 7 . 5 ] 2 -7 H9O. (30131] 

Many other frustrated spin models in one, two or three 

In the present paper we combine the ideas of Shas- 

dimensions are known which have also dimer-singlet prod- 
try and Sutherland [12] and our recent findings on ex- 

uct states as groundstates, see, e.g., [18I19I20I21I22] . A 

act trimerized singlet product ground states (TSPGS's) 

systematic investigation of systems with dimerized eigen- 

for ID integer-spin Heisenberg systems |32] and discuss 

states can be found in [23] ■ Note that these dimer-singlet 

such TSPGS's on a two-dimensional modified Shastry- 

product groundstates have gapped magnetic excitations 

Sutherland square-lattice model. Section [2] shortly reca- 

and lead therefore to a plateau in the magnetization at 

pitulates the theory of TSPGS's and section [3] defines the 

m = 0. Recently it has been demonstrated for the ID 

modified Shastry-Sutherland model and its finite realiza- 

counterpart of the Shastry-Sutherland model [19121122] . 

tions that will be analyzed in what follows. We have con- 
that more general product eigenstates containing chain 

centrated in our numerical studies on the size of the gap 
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for the exact ground state for finite lattices of N = 12 (for 
spin quantum numbers s = 1, s = 2), as well as N = 18 
and N = 24 (for s = 1) and on the magnetization curves 
for selected values of J 2 , see section l4~T1 The analytical re- 
sults in section l4~2n mainlv concern upper and lower bounds 
of the gap function. These results depend on a slightly 
generalized statement and proof of the gap theorem, first 
formulated in |32| . which is done in appendix [21 Finally, 
appendix [B] contains exact results on classical magnetiza- 
tion curves for the model under consideration. 

2 Exact ground states 

The anti-ferromagnetic uniform spin trimcr 

Hi = J(s ■ si + s ■ s 2 + si ■ s 2 ) (1) 




Fig. 1. The modified Shastry-Sutherland model on the 
decorated square lattice for TV = 24 sites (periodic condi- 
tions imposed) used for exact diagonalization. 



has, for J > and integer s, a unique S = ground state, 
denoted by [0, 1,2], with ground state energy 



--Js(s + 1) 



The corresponding product state 

M 



£0,11, £2] 



will be an eigenstate of a system of Af coupled spin trimers 
indexed by i = 1, . . . , Af with Hamiltonian 



for all 1 < i < j < Af and 5, e = 0, 1, 2, see [35]. Moreover, 

(2) fl3j will be a ground state of Q, a TSPGS, if the intra- 
trimcr coupling is almost uniform and the inter-trimer 
coupling is not too strong [32] . The domain of coupling 

(3) constants where this is the case will be called the "TSPGS- 
rcgion" . 



iejS 



(4) 



if and only if the coupling between different trimers is 
"balanced" in the following sense: 



J, 



i5,j5 



If the system of trimers has a periodic lattice structure, 
the difference AE between the energy of the first excited 
state and that of the ground state can be shown [35] to 
be bounded from below independently of the system size. 
(5) In other words, the TSPGS is "gapped". 




Fig. 2. Two finite decorated square lattices of A" = 12 
and N = 18 sites used for exact diagonalization. 

3 The model 

We consider the inter-spin Heiscnbcrg Hamiltonian on a 
decorated square-lattice (see figured]). It results from the 
well-known Shastry-Suthcrland model by replacing its di- 
agonals by equilateral triangles with uniform intra-trimer 
interaction strength J\ > 0. The set of triangles is di- 
vided in a bi-partitc fashion into two disjoint subsets of 
triangles of type I and type II, corresponding to diagonals 
with positive slope resp. negative ones, see figure [T] Each 
triangle of, say, type I is surrounded by four triangles of 
type II and connected to each of them with three bonds 
of strength J 2 . 

It follows that the inter-trimcr coupling satisfies the bal- 
ance condition ((5]) and hence the theory of TSPGS's ap- 
plies. In particular, two questions arise which will be ad- 
dressed in the following sections: What is the size of the 
TSPGS-region and of what kind are the lowest excita- 
tions? The latter question is also connected to the issue 
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of magnetization plateaus which will be shortly discussed 
below. 

4 Results 

4.1 Numerical results 

In what follows we set J\ = 1 and consider J 2 as the 
variable bond strength. To study the region where the 
TSPGS is the ground state of the model ^ we use the 
Lanczos exact diagonalization (ED) technique. Since for 
spin quantum numbers s > 1/2 considered here the size 
of the Hamiltonian matrix grows much faster with system 
size N than for s = 1/2, we are restricted to finite lattices 
of N = 12, 18 and 24 for s = 1 and N = 12 for s = 2. The 
largest lattice is shown in figure [TJ whereas the smaller 
lattices are shown in figure [5J Although the criterion for 
the existence of TSPGS's (see section \S§ arc fulfilled, we 
have to mention that for the small lattices of N = 12 and 
N = 18 the exchange pattern of the J\ diagonal bonds in 
the squares do not match to the infinite system. Neverthe- 
less, we have included the data for N = 12 and 18 to get 
an impression on finite-size effects and on the influence of 
the spin quantum number s. 

According to [32] the TSPGS is gapped. Hence we use 
the spin gap, see figure [3] to detect the critical points Jf 1 
and Jf 2 , where the TSPGS gives way for other ground 
states. We find for s = 1 the values J 2 cl = -0.570, -0.578, 
and -0.587 and J 2 cl = 0.434,0.446, and 0.454 for N = 
12,18, and 24, respectively (cf. figure j3ja)). For s = 2 
and N = 12 we have Jf 1 = -0.400 and Jf = 0.322, 
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cf. figure EJb). These values lie between the upper and 
lower bounds which will be derived for J^ 1 an d J% 2 m the 
next section for N — > oo. The nature of the lowest ex- 
cited state depends on J^. Around J 2 — it is a triplet 
state with strong antiferromagnetic correlations along the 
trimer bonds and weak correlations between the trimers. 
Near Jf 1 the lowest excitation is a fcrrimagnetic state, 
i.e. the total spin is S = As/3 and the system splits 
into two ferromagnctically correlated sublatticcs contain- 
ing on the one hand the 2A/3 square-lattice sites (i.e. 
sites 0, f , . . . , 15 in figure [T]) and on the other hand the 
A/3 additional sites (i.e. sites 16, 17, . . . , 23 in figure [T]). 
The spin correlations between both sublattices are anti- 
ferromagnetic. The ferrimagnetic state is the ground state 
for —1.5 < J2 < J2 1 ■ Near J| 2 the lowest excitation is 
a collective singlet state with strong correlations along 
all bonds, and, this state becomes the ground state at 
•h = J 2 ''■ 

It is well known that the magnetization curve of the 
Shastry-Sutherland model (as well as that of the corre- 
sponding material SrCu2(B03)2) possesses a series of pla- 
teaus, see, e.g., [27 33 34 35]. Motivated by this, we study 
now briefly the magnetization curve M(h) (where M is 
the total magnetization and h is the strength of the ex- 
ternal magnetic field) for the considered model for s = 1 
using ED for TV = 18 and A = 24 sites. ED results for the 
relative magnetization m = M/M sat versus magnetic field 
h for N = 18 and A = 24 sites are shown in figure H^i. 
Again the finite-size effects seem to be small. Trivially, in 
the limit J2 — the m(h) curve consists of three cquidis- 



LU 



LU 

< 




-0.6 -0.4 -0.2 0.2 0.4 0.6 

Fig. 3. Numerical exact data for N = 12, 18, and 24 
(symbols) as well as upper (black solid line) and lower 
bounds (red solid line) for the excitation gap AE. (a) spin 
quantum number s = 1; (b) spin quantum number s = 2. 
Note that the labels S = 1, S = 0, S = 2 A/3 (fcrri), 
and S = 8 (ferri) characterize the total spin of the excited 
state. 
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m=1/3, N=18 

m=2/3, N=18 

1 1 ^ 1 1 1 1 

-0.6 -0.4 -0.2 0.2 0.4 

J 2 

Fig. 4. (a) Magnetization curve m(h) for selected values 
of J 2 and s = 1 (thick lines N = 24, thin lines N = 18); 
(b) Plateau widths Ah of the m = 1/3 and the m = 2/3 
plateaus as a function of J2 for N = 24 and TV = 18 and 
s = 1. 
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tant plateaus and jumps according to the magnetization 
curve of an individual triangle. Switching on a ferromag- 
netic inter-triangle bond J2 < the general shape of the 
magnetization curve is preserved. However, the saturation 
field as well as the end points of the plateaus decrease al- 
most lineraly with J2 and become zero at J2 = 1.5, where 
the ground state becomes the fully polarized ferromag- 
netic state. 

In case of a moderate antiferromagnetic inter-triangle 
bond J 2 > the plateaus at m = 1/3 and m = 2/3 still ex- 
ist, however the discontinuous transition between plateaus 
becomes smooth. Note that a m = 1/3 plateau was also 
found for the standard Shastry-Sutherland model |34l35j . 
The plateau widths Ah of the m = 1/3 and m = 2/3 
plateaus in dependence on J2 is shown in figure 0Jd. Ob- 
viously, both widths shrink monotonously with increasing 
of I J2 1 . If J2 approaches the critical value J^ 1 we find in- 
dications for additional plateaus, e.g., at m = 5/6. Note, 
however, that our finite-size analysis of the plateaus natu- 
rally could miss other plateaus present in infinite systems, 
see, e.g., the discussion of the ED data of the m(h) curve 
of the standard Shastry-Sutherland model in |36j . Hence, 
the study of the magnetization process of the considered 
quantum spin model needs further attention based on al- 
ternative methods. 

One might expect that the presence of these plateaus 
and jumps may be linked purely to quantum effects be- 
cause they are often not observed in equivalent classical 
models at T = |25I37I38I39] , However, for the present 
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4.2 Analytical results 



7 



4.2.1 s = 1 




Fig. 5. Two possible subsystems of the modified Shastry- 
Sutherland lattice, see e. g. figure Q] The upper one, Hq, 
consists of two coupled triangles; the lower one, Hg of 
three triangles. 



model the plateau at m = 1/3 survives in the classical 
limit for J 2 < as we will show in appendix [Bj 



In order to obtain analytical results about the TSPGS- 
region we have adapted theorem 3 of [35] to the present 
situation. A slightly more general version of this theorem 
is stated and proven in appendix [21 It yields lower bounds 
for the gap AE of the form AE > /(4J) and the TSPGS- 
region in terms of properties of simpler spin systems of 
which the lattice can be composed, see figure [3] These 
subsystems are chosen here as systems isomorphic to Hq, 
see figure [SJ consisting of two neighboring triangles. For 
s = 1 the gap function x = SqE = /(J), J = of 
Hq is obtained as a special case of equation ([7]) given be- 
low. This yields the corresponding bounds for the TSPGS- 
region ( J cl , J c2 ) 

jd 3_V73 _ _ Q 3465 1 jc2 fo s = 1 /gj 
16 4 w 

The function SqE = f(J) according to ([7]) also provides 
an upper bound for the gap function of the lattice, since it 
represents the energy of a state orthogonal to the TSPGS, 
albeit not an eigenstate of H . This bound is very close to 
the numerically determined gap function in the case of 
N = 12, see figure [3^,, but considerably deviates in the 
cases of N = 18 and TV = 24. This indicates that, in gen- 
eral, the lowest excitations of the lattice arc different from 
the excitations of Hq. 
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4.2.2 General s 

It is possible to analytically calculate the energy of the 
lowest excitations of H§ for general integer s. The corre- 
sponding gap 5qE = x = /(J) is obtained as the lowest 
root of the following cubic equation 

- (x - 4) (re - 2) (a; -l)-(x- l)(2x - 5) J 

+ (1 - 3r- x + rx)J 2 +rJ 3 = (7) 



<5 6 (E) 
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where we have set r = s(s+l). From this result one derives Fi S- 6 - Lower bounds of thc scalcd S a P function 
the lower bound 



AE > /(4J) for general s 



(8) 



6eE(j), j = ys(s + 1) J of the modified Shastry- 
Sutherland spin lattice for s = 1,...,10 (thin curves) 
obtained from Eq. (Q- The curves approach the asymp- 
and a lattice of arbitrary size, sec theorem 1 in appendix^] totic (HH]) for s -> oo (thick red curve) which has a simple 
adapted to the system under consideration. The corre- quadratic approximation dTTJ> (thick green curve), 
sponding curves are shrinking in J-direction with increas- <5o(E) 
ing s and yield inner bounds for the TSPGS-rcgion ( J cl , J c2 ) 
of thc form 



j c1 < < o < j£ 2) < j c2 , 



(9) 



see figure [5] (green curves). Upon scaling w. r. t. the new 
variable j = y/rJ the graphs of ([7]) asymptotically ap- 
proach the curve given by 

;2 (a;-4)(i-2)(a:-l) 



M 08 




/ml 06 




ft 1 0.4 




i/m I 0.2 





J 



16(x- 3) 



with Taylor expansion 



x = 1 



32 



see figure[ni Hence J 1 / assumes for s — > oo asymptotically 



(10) Fig. 7. Upper bounds of the scaled gap function 
SoE(j), j = \J s(s + 1) J of thc modified Shastry- 
Sutherland spin lattice for s = 1,...,10 (thin curves) 

(11) obtained from Eq. (|15l) . Thc curves approach the asymp- 
totic p7[) for s — > oo (thick red curve) which has a simple 



the form 



quadratic approximation (|18[) (thick green curve) 



\J 



(0 



1,2 



(12) 
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(J cl , J c2 ) of the form 



In order to obtain close upper bounds g(J) of the gap 
A(E) in the case N > 18 we calculate the energy of a 
certain (degenerate) state that involves three triangles for 
arbitrary integer s, say, one triangle of type I and two 



4 X) < j c1 < o < r 2 < , (16) 

see figurc[5](red curves). Upon scaling w. r. t. the new vari- 
able j = y/rJ the graphs of (fl~5j) asymptotically approach 
the curve given by 



neighboring triangles of type II, see figure This state j 2 = 2±E 221^. rZl^ il 

J 7a; - 15 ' v ' 



3(ar - 3)(x-2)(x - 1) 
with Taylor expansion 



is obtained as an exact eigenstate of Hq, which is the full 
Hamiltonian H , restricted to a 4 3 = 64— dimensional sub- 
space spanned by product states of the form 

x=l--f + 0(f), (18) 



b k , i,j,k = 0, ...,3. (13) 



see figure [3 Hence Jfj assumes for s — > oo asymptotically 
the form 



The 4> n live in the (2s + 1) -dimensional Hilbcrt spaces be- | t(*)i / ° -_i o (in\ 

1 u '~y 5s( 8 +iy l - l > 1 - ^ 

longing to one of the three triangles. <po = [0,1,2] denotes 
the TSPGS of the corresponding triangle and 

(«), 

, i — 1,2,3 , (14) Although these curves constitute only upper bounds 



I £0 9o 



of the true gap functions, the comparison with the nu- 

(i) 

where sX is the i-th component of the spin operator Sn ,cat 10 i »t ^ ^ 

~ u ~ merical results for TV = 18 and TV = 24 reveals a close 

pertaining to the spin site number 0, an arbitrarily chosen . . , , . „ „ . 

approximation to both curves, see figure [3| Ihis supports 

spin site of the corresponding triangle. The gap function 



of Ho will be denoted by x = SqE = g(J) and constitutes 
an upper bound for A(E). It has the following implicit 
form, using r = s(s + 1): 



our conjecture that (JTSJ) indeed may serve as an analytical 
approximation of the gap functions for large N and arbi- 
trary integer s. This would mean that the excitations from 
the TSPGS can be viewed as local excitations essentially 

= -12(x - 3) 2 (x - 2){x - 1) - 6J(x - 3)^ - l)(4x - 9) conccntratcd 011 thrcc neighboring triangles. Numerically 

determined spin correlation functions seem to be in accor- 

+J 2 {x - 3) (9 - 9a; + 4r(7a; - 15)) + l6J 3 r(2x - 5) . 

dance with this conjecture. Of course, the corresponding 

(15) 

excited state will be largely degenerate due to the trans- 
Again, the function g belongs to the lowest branch of (fT5"j) . lational symmetry of the lattice. We expect an almost flat 
The corresponding curves are shrinking in J-direction with k-dependance of the energy band E(k). This expectation 
increasing s and yield outer bounds for the TSPGS-rcgion is also supported by our numerical results. We have found 
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Fig. 8. Exact bounds for the TSPGS-region ( J cl , J c2 ) 
for s = 1, . . . , 10 of the form < J cl < < < j[ 2) < 
J c2 < J^ 2 ' . These are derived from Q (green curves, inner 
bounds) and (|15[) (red curves, outer bounds). In the classical 
limit s — > oo the TSPGS-region shrinks to zero according to 
(fT5)) and (HU). 

that the lowest excitations close to J = have the to- 
tal spin quantum number S = 1 in accordance with our 
model. 

In the case iV = 12 where we have performed numerical 
calculations for s = 1 and s = 2 it is not possible to put 
a subsystem of type Hg into the lattice and the above 
results do not apply. However, an analogous method can 
be applied to two coupled triangles of type Hq and yields 
an upper bound of the gap function of the form 

AE < ^(18 - 3J - \/9(J - 2) 2 + 96J 2 s(s + 1)). (20) 

The numerically determined gap together with the bounds 
(|8|) and (|20| is represented in figure [3] b. 
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A Proof of the gap theorem 

In order to prove the existence of an energy gap between 
the TSPGS and the first excited state we will adapt the 
analogous proof given in 32! to the modified Shastry- 
Sutherland model considered in this article. The gap the- 
orem will be formulated in a slightly more general frame- 
work. 

For any positive integer £ let Z^ denote the set of in- 
tegers modulo £, such that n + £ = n, and 

Z L = Z Ll x . . . x Z Ld (21) 

a standard d-dimensional lattice with total size 

(/ 

Af = Y[L £ . (22) 

e=i 

We consider an index set C on which the additive group 
Zl operates effectively, i. e. without fixed points except for 
the neutral element. Let c be the number of correspond- 
ing equivalence classes (orbits), c = \£/Zj_,\. Each orbit is 
isomorphic to Zl ; if we select an index from each orbit we 
obtain a bijection ip : {1, . . . , c} x Zl — > C which we call 
a "(global) trivialization" in analogy with the correspond- 
ing term in the theory of fibre bundles. 
We will fix t different trivializations 

<p v : {1, . . . , c} x Z L — ► C, v = l,...,r. (23) 
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Fig. 9. The Hamiltonian of the (modified) Shastry- 
Sutherland spin lattice can be written as the sum of four 
classes of simpler Hamiltonians of the kind Hq , sec figure 
[3 which are indicated by the colors black, red, green, blue. 
The subsystems of each color are isomorphic via trans- 
lations T of the lattice. They can be transformed into 
each other by 90° rotations R about the mid-points of the 
empty squares. Note that there are two types of diagonals 
(trimers), one with positive slope (type I) and the other 
with negative one (type II). Each diagonal belongs to four 
subsystems of different color, hence the coupling constant 
Ji of the subsystems, see figure has to be divided by 
the factor 4 in order to obtain a total Hamiltonian with 
the coupling constant J\. 

W. r. t. these trivializations the total Hilbcrt space can be 
written as a tensor product space in the following form: 

n = ^n l = (g) n vAe , k) (24) 

ieC (£,k)g{l,..,c}xZ L 
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= (g) (g) «M*,k) (25) 
kez L e=i,...,c 

= (g H v>k , u = l,...,r. (26) 

k(EZ L 

We will also write 

W„,k = (g) Hi where S uM = {<^(£,k)|l < I < c} . 

(27) 

We will explain these definitions in the case of the 
modified Shastry-Sutherland model in the finite realiza- 
tion of figure [TJ It is a d = 2-dimensional lattice with 
Li = L 2 = 2, hence N = L ± L 2 = 4. The 8 different 
triangles of figure [T] correspond to the indices i G C, 
hence C can be viewed as the set of triangles and be 
identified with the set of numbers of their upper cor- 
ners {16, 17, . . . , 23}. Note that C is not isomorphic to 
the underlying spin lattice which has 24 sites. The Hilbert 
spaces Hi, i € C happen to be isomorphic and of the 
same dimension (2s + l) 3 . In general, it is not necessary 
that all Hi are isomorphic; e. g. we could have spin sys- 
tems composed of trimers and dimers. Zl operates on C 
in a natural way by means of translations. Type I tri- 
angles (18, 19, 22, 23) cannot be transformed into type II 
triangles (16,17,20,21) by means of translations. Hence 
we have, in this case, c = |£/Zl| = 2. In figure [5] two ad- 
jacent triangles of different types are coupled together in 
order to form sub-Hamiltonians of the kind Hq summing 
over all possible translations of it. These sub-Hamiltonians 
are in a 1 : 1 manner characterized by trivializations 
of the kind considered above: We denote the two adja- 
cent triangles we started with by (1,0) and (2,0), and 
all other pairs which are translations of these will be de- 
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noted by (1, k) and (2, k) where k runs through Zl- Since 
the set £ is exhausted by this construction we have ob- 
tained a trivialization in the sense of (|23p . i. e. a bijection 
p„:{l,2}xZ L — >£,v = 1,...,t. In figure [9]t = 4 dif- 
ferent trivializations are shown and indicated by different 
colors. 

Returning to the general case we will identify all fac- 
tor spaces belonging to subsystems of the same type by 
means of a certain product basis {I7)}, 7 : £ — > N in % 
such that i G £ 1 — y G {0, . . . , dim;}, where dim^ de- 
notes the dimension of Hi- Correspondingly, the unitary 
translation operators T m , m G Zl are defined by suitable 
permutations of the product basis. W. r. t. a trivialization 
<p v this definition assumes the form 

(T m 7 )(M4 n )) = ibPvV, m + n)) . (28) 

The T m , m € Zl form the abelian translation group T 
whose characters arc of the well-known form 

x(T m ) = cxp(2^m-k/A0, (29) 
where m = (mi, . . . , m<j), k = (fci, . . . , fc<j) G Zl(30) 

The total Hamiltonian H is assumed to be a sum of sub- 
Hamiltonians of the form 

r r 

ff-EE ^EE T^HvT^ , (31) 
f=ikez L y=ikeZ L 

where -ff„.k is defined on the "supporting factor space" 

"H^k, see (|26[) . and extended as the identity operator on 

the remaining factor spaces to the total space H. Of course, 

it suffices to postulate this only for H u — H v q. 
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As a consequence of (|3"Tj) we note that total Hamilto- 
nian will commute with all translations: 

[H, T m ] = for all T m G T . (32) 

Moreover, we assume the following: 

Assumption 1 For all i G £ let <Pi G Hi be normalized 
states such that 

<P = ®*i (33) 

is a ground state of H u for all v = 1 , . . . , r which is unique 
on the factor space H. Vi q. We set 

H u $ = El Q) $ (34) 

and denote by 

E<V = E^ + 8 U , 5 V > (35) 

the next-lowest energy eigenvalue of H v . Moreover, <& is 
assumed to be invariant under translations, 

T m & = $ for all T m eT . (36) 

Then the gap theorem can be formulated as follows. 

Theorem 1 Under the preceding definitions and assump- 
tions <P will be the unique ground state of H with eigen- 
value Eq = Af Y^1=i and the next-lowest eigenvalue 
of H satisfies E\ > Eq + Y^l=i &v 

The existence of a gap follows since X^=i &v > i s m ~ 
dependent of the size M of the lattice. In the special but 
important case where all H v , v = l,...,r are unitar- 
ily equivalent, we write E' 1 ' = E^ + S and conclude 
E a =NtE^ and E t > E a +tS. 
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Proof of theorem [JJ The first claim (except uniqueness) = E^ (a\ W^' k \a) (42) 

follows immediately from assumption [1] and the fact that, / \ 

J2 4 0) (0|^' k |0) + Y,Ei a) {a\W^\a) (43) 

I/,k V a>0 / 

> J2 ( 4 0) (o|^' k |o) + (4°) + y^H< 



\a) 

a>0 



due to [361 ^ is also a ground state of all i?„.k with the 

same eigenvalue E { °\ k e Z L and A/" £° =1 £^ 0) being an > \^ ( r(o) / ft lW/f' k ln\ 4- CrW + ^..^ \^ ln\W^ 
obvious lower bound of H. 

Let £ % be the eigenvector of H belonging to the next- 
lowest eigenvalues E\ > Eq. We first note that W _L <P 
follows in the case E\ > Eq and can be arranged in the (0|^V 10) < 1 - • (45) 



(44) 



Lemma 1 



case Ei = Eq (which we cannot exclude from the outset) 
by choice of W. Moreover, due to (|3"2"]l we may choose \P to 

Proof of lemma [TJ It suffices to consider the case k = 

be a common eigenvector of all translations, 

0. We again consider the product basis {I7)} introduced 
^ = x(T m ) & for all T m g T , (37) above and write the c quantum numbers ■y(ip u (l 7 0)), I = 

1, . . . , c, at the first c places of the string I7) = 

where x(^m) is of the form (|29l) . 

I rii, ri2, ri c , ...). It follows that the ground state in 

Our aim is to show i?i > £0 + X^=i Let I 01 ); a = 

is denoted by a ket |0, 0, . . . , 0) consisting of c zeroes. 

0, 1, 2, . . . denote the eigenbasis of H v y_ in "H„ ; k- Further we 

We conclude 

arrange the eigenbasis such that |0) = ® ieS k holds, 

see (EZJ). The corresponding eigenvalues of H uM are de- (0, . . . , 0| W£'°|0, . . . , 0) = V] |(f |0, . . . , 0, X) | 2 = s 
noted by Ev*' 1 , in accordance to the notation E^ and E^ (46) 
introduced above. \a, K) denotes a corresponding product anc | 
basis in %, where K stands for some multi- index of quan- 



tum numbers. Moreover, we consider the reduced density 



1 = TiW£° (47) 



operator^ in ^.defined by the partial trace £ l« ■ • • , 0, n c+1 , n c+2 , . . .)|* (48) 



ra c+ i,n c+ 2, 



(a\W^\(3) = *£{a,K\V){V\p,K) . (38) 



+ ]T K!?|m,n a ,...)| 2 (49) 



ni,n 2 ,... 



^ = S + Si ■ (50) 

Then we conclude 

The first sum so in (|4"T|) runs through all sequences 
E x = (V\H\V) = ^2(W\H Utk \W) (39) 0, . . . , 0, n c+1 , n c+2 , . . . excluding the value n c+1 = n c+2 = 

i/,k 

_ Tj^y rj^. ^ W^'^j (40) • ■ • = ^' smce (^1^) = 0- Equivalently, we will say that it 

runs through all states ip = |0, . . . , 0, n c +i, n c +2, • • •) S 

= £ Tr(4«)|a)(aK' k ) (41) 

„ k v y The second sum s\ in (|48[) runs through all sequences 
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77,1,77,2,... except those with n% = 77,2 = ■ . ■ = n c = 0, ^ Aeyl t> = is impossible since it would imply that t> = 

or, equivalently, through all states ?/> = |ni, n%, . . .) G B\. for all A € A and hence 1 = sq + Si = 0. 

Thus the total sum in (|47I48[) runs through an orthonor- From (|55[) wc infer 

mal basis £ = Bo U #i of W = {tp € H\(tp\^) = 0}. -, -, 

We consider on B the equivalence relation ipi ~ ip 2 ^ S ° S ° S ° 

V'l = T m ip2 for some T m € T, and denote by A — B/^ and 

the corresponding set of equivalence classes or "orbits" . Mil 
Due to (|37[) all states -0 in the same orbit A yield the same 7V • /v 

value 51 ^ ' (58) 

which concludes the proof of the lemma. □ 

t x ^\m)\ 2 = mT m ^)\\T m ^r . (si) 

For each orbit A G A let N x = |A| denote its length. For To complete the proof of theorem[T]we use again the eigen- 
most orbits we have N x = TV, but in general N x will basis of H vM and write 
be a divisor of N '. For example, if d = 1, L\ = N = Q 



and |1,2,3,1,2,3) G A then N x = 3. We define N { x k) 
|A n Bk\, k = 0, 1, and obtain the following equations: 



A^a = ^i 0) + N ( x 1] , (52) 



so = x>M 0) . ( 53 ) 

AS/1 

Sl = ^t A ivf ) . (54) 



1 = Tr = (0|W^ k |0) + £ («|W£ k |a) = s + Sl . 

a=l,2.... 

(59) 

Then wc rewrite (jUJ) as 



Si >^(4 0) So + (4 0) +^)si) (60) 

= X)(4 0) +^*i) (61) 

i/,k 

= lo+^X]^ s i >-Bo+^viX]^ ( 62 ) 

T 

= £ +^<5„, ( 63 ) 



AeA 

Let ^ = |ni, ?i2, 713,714, .. .) G A. Note that at least one 
n ji 1 < j < cj\f must be non-zero since ip ^ # = 
|0, 0, . . . , 0). Hence at least one translation of V belongs whcre we haye ugcd ^ which is cquivalent to i emma [T] 
to B\ , namely that where j is shifted to one of the first c ^ 
places. To show this in detail wc write rij — j(ip u (£, k)), 1 < 

C, k G Zl- It follows that T-k^ = |TOi , 7Tl2, TO3, 7714, ...) G t J l-i-l ,i . ,i j-n j 

— ' u Kri-L)-!)aj4)/ j n orc i cr to apply the gap theorem to the modified 

A and mi = n,-, ^ 0. Thus > 1 and hence < C i ± c j-u i j i j. t n • r , i 

a — a — Shastry-Sutherland lattice we write lor the energies of the 

N X - 1 < AT - 1 which for E A£ a h > implies subsystems H V)k 

T Q -Af—[- (55) ^ 1 )=^°)+5(|,J 2 ). (64) 
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This has the consequence that the total Hamiltonian will 
correspond to the coupling constants J\ = 1 and J2 since 
each triangle is contained in four different subsystems, 
see figure [9] Since 5 is a homogeneous function of J\ , J2, 
i. e. <5(aJi,aJ 2 ) = aS(Ji, J2), we may write 6(j, J2) = 
j 6(1, 4J 2 ). Hence the gap theorem implies (note that r = 

E > £i+4i(5(l,4J2) = £i+<5(l,4J 2 ) . (65) 

Thus the lower bound of the gap is simply obtained by 
shrinking the graph of the gap function 6(1, J 2) of Hq into 
^-direction by a factor 4. 

B Classical ground states 

It follows from the general theory [35] as well as from 
our special results (jTU)) and (JT7J) that the classical modi- 
fied Shastry-Sutherland model possesses no TSPGS's ex- 
cept for J2 = 0. Nevertheless it is possible to analytically 
obtain the classical ground states for J2 7^ and arbi- 
trary magnetic field h and from these the magnetization 
curves. Typically in the classical limit s — ¥ 00 the magne- 
tization curves at T = are smooth and do not exhibit 
plateaus or jumps [25I37I38I351 . An exception to this rule 
is, e. g., reported in [3D]. Hence it is remarkable that the 
classical modified Shastry-Sutherland model possesses a 
plateau at a magnetization of AI = ^Af sat and a jump at 
h = 0, J2 = —1/2, as will be shown in the sequel. 

For sake of simplicity we assume a quadratic square 
lattice of L x L squares, where L > 4 is some multiple of 
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Fig. 10. Magnetization curves for the classical modified 
Shastry-Sutherland lattice for different values of J = — J 2 
between and 0.6. The magnetization plateau at M = 
■|Af sat (indicated by red color) separates states of phase 
II (blue curves) from those of phase III (black curves). 
Full magnetization M = Af sat is reached in the region 
indicated by green color. The phase boundaries can be 
analytically calculated, see (j72j) - (f76|) . 

4. It hence contains I = ^- triangles fi = 1, . . . , £. Let 
s ai,0j s m,1i s m-2 denote the three (unit) spin vectors corre- 
sponding to such that s^o corresponds to the "out-of- 
plane" spin site, see figure [JJ and 

S A , = s M , + s^i + s M;2 , [-1 = 1,..., I (66) 

denote its total spin. is uniformly coupled to two ad- 
jacent sites which belong to neighboring triangles with 
strength J2. As usual, we write the Zeeman term in the 
Hamiltonian as —hM = —h (X^i=i where h is the 
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Fig. 11. A typical classical ground state of phase II ac- 
cording to ([oT)) - (fT0)) . The numbers /i which are attached 
to the spin vectors s M refer to the spin sites of figure [TJ All 
spin vectors of the in-plane spins lie on a cone with open- 
ing angle 6, whereas the out-of-plane spins have a polar 
angle 

strength of the (dimcnsionslcss) magnetic field. We will 
confine ourselves to the ferromagnetic case J2 < 0, which 
shows the most interesting features. In the AF case J2 > 
the magnetization curves are almost linear until they reach 
the saturation domain. 

As one can see in figure \W\ there are, besides the fully 
aligned state with M = M sat = 3£, exactly three different 



phases. Phase I which forms the magnetization plateau at 
M = ^Afgat = £ is given by the uud-st&te, i. e. in each 
triangle the two in-plane spins point into the direction 
of the magnetic field ( "up" ) and the off-plane spin in the 
opposite direction ("down"). 

The two other phases II and III have spin vectors of the 
form 



^sin6> cost^ ^ 



sine/ siny^, 
cos# 



and 



1,2 (67) 



V 



sin $ cos tp^ 
s (i0 = sin sin^o 1 fi = l,...,£. 
cos$ 

Ground states of phase II live in the domain < M < 
^M sat and are visualized in figure 111! Their azimuthal 
angles ip^ j , i = 1,2 assume 8 different values 



Vt^ : = n- ± 4>, n = 0, . . . ,3 



which depend on a parameter 4>, whereas 



^0 = n—, 7i = 0, ... ,3 . 



(69) 



(70) 



Upon a translation from, say, the triangle (0, 7, 16) to 
(4, 9, 19), see figure [1] the azimuthal angles (|6"9")l and (fTOf 
are shifted by an amount of — | . Hence the states of phase 
II are characterized by a wave number of k = -| (the minus 
sign does not matter). If the magnetic field h approaches 
the left hand boundary of the M — ^M sat plateau, the 
states of phase II become the uud-state. Phase III is con- 
fined to a magnetization M satisfying ^M sat < M < M sat . 
The corresponding states have a wave number k — n, since 
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their azimuthal angles satisfy tp^i e {0, n} for (j, = 1, . . . ,£ Actually, the phase boundaries, displayed in figure [TU] 
and i = 1,2,3. More precisely, the spins in figure [T] with by thick lines, can be given in closed form. We will write 

T 2 

the numbers 0, 2, 5, 7, 8, 10, 13, 15; 18, 19, 22, 23 have an az- J = — J2 and remind the reader of Ji = 1 and I = -=5- = 

imuthal angle of ip — 0, and the remaining spins of ip = ir. |M sat . The curve for h — > (thick blue curve in figure ITUf 

It is obvious how to generalize the states of phase I, II, consists of three parts, the phase II part 
III to infinite modified Shastry-Sutherland lattices by pe- 
riodic continuation. The (semi-)analytical treatment of 

these states can be based on the ground state equation, the P nase M P ar t 



M n (J) = IJ V1 + 2J, < J < - , (72) 



see 







Mm(J) = ^ V J(l + 2J)(3 + 2J), -<J<-, (73) 

and a jump at J = | from Mi = to M2 = L 
fj, = l,...,N. (71) The saturation field h sa t(J) (thick green curve in figure 
[TU|) is given by 



nt wv, u r v ^ 1 t ^sat(J) = i (3 - 11J + ^3(3 + 10J + 19J 2 )) , < J < ^ 

Mote that here Greek indices /i, u = 1, . . . , N do not num- 2 V 12 



ber triangles but all spins of a finite spin lattices. The 
ft M denote Lagrange parameters due to the constraints 
s p, ' s /j . = 1- They can be eliminated by writing j M x s M = 
0, /i = 1, . . . , N. These equations, together with the cor- 
responding periodicity properties of phase II or III states, 
are sufficient to determine the unknowns 9, •& and, for 
phase II states, 0, as solutions of certain algebraic equa- 
tions that contain J2 and h as parameters. We have cal- 

. ._. Finally, the plateau at M = \M sa , t (bounded by thick 

culated tne magnetization curves of figure [lUj by means 



(74) 

The value /i sa t(0) = 3 is part of the linear magnetization 
curve for J2 = (thick black line in figure [TO]) M(h) = £h 
which corresponds to the magnetization of a uniform AF 
triangle. 

The two boundaries (f73"| and (|74p meet at the point h = 

0, J = §, M = M 8at = 3*. 



red curves in figure [T0"|) is given by the inequalities hjj(J) < 
h < hi 1 1 (J) where ha (J) is the lower positive root of 



of numerical solutions of these equations and checked the 
results by a direct numerical calculation of the ground 
states for different h and J. For h = we found a large /i(/i— l) 2 +4(/i— l)(2/i-l) J+(19/i— 16) J 2 +16 J 3 = (75) 



in the interval < J < i and 



degenerate set of ground states with total spin S varying 
between and some >!5 max . In the limit h — > only the 
states with S = S max are obtained as ground states with hju(J) — — ^1 — 9 J + y/l + 14 J + 17J 2 ^j , < J < — . 
finite magnetization M = S'max. (76) 
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Obviously, the graphs of hjj and hju intersect at the two 
points J = 0, h = 1 and J = |, h = 0. 
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